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AMENDMENTS T Q THE CLAIMS 

1. (Original) A machine-processing method for computing a property of a 
mathematically modelled physical system, the steps comprising: 

a) reading, via a machine processing unit, input data including a value for each 
identified ordered coefficient of a first polynomial p(x) representing said property, said 
polynomial p(x) being expressed asp(x) - I^Pj • x*) where j=0 to n, a value of a quantity 
x, a value of a predetermined Xi, and a value of a predetermined p(xj) correspondingly 
paired with said predetermined x*; 

b) building, via said machine processing unit, a value of a second polynomial c(x) 
having ordered coefficients, said second polynomial c(x) being expressible as: c(x) - 
2(C k • x*) where k=0 to (n-1) so that said first polynomial p(x) is expressible as: 
p(x)=p(xO+{x- Xi } • c(x), comprising the steps of: 

i) determining, via said machine processing unit, a value for each ordered 
coefficient of said second polynomial c(x) by generating a plurality of machine 
processing unit signals to determine each said ordered coefficient of said second 
polynomial c(x) from: C k - KJVnj> • x";) where j=0 to (n-14c); 

ii) determining, via said machine processing unit, a value of said second 
polynomial c(x) by generating a plurality of machine processing unit signals to 
determine: c(x) = 2(C k • x*) where k=0 to (n- 1); 
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c) constructing, via said machine processing unit, a value of said first polynomial 
p(x) by generating a plurality of machine processing unit signals to determine: 
p(x)=p(Xi)+{x-Xi} • c(x); and 

d) ourputting, via said machine-processing unit, said value of the first polynomial 
p(x) representing said property of the mathematically modelled physical system. 

2. (Original) The machine-implementable method of claim 1, wherein a 
difference between x and x ; is sufficiently small to achieve a desired accuracy of a final 
computation of said machine representation of a numerical value of said first polynomial 

3. (Original) The machine-implementable method of claim 2 wherein the step of 
reading said input data comprises reading, via said machine processing unit, said input 
data from a machine-readable medium. 

4. (Original) The machine-implementable method of claim 3 wherein said 
ordered coefficients of said second polynomial c(x) are computed from a mathematical 
expression derivable from: C k = ZtP^i+j, • x j i ) where j=0 to (n-l-k). 

5. (Original) The machine-implementable method of claim 4 wherein said 
mathematical expression is a mathematical recurrence expression. 
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6. (Original) The machine-implementable method of claim 5 wherein said 
mathematical recurrence expression is a forward mathematical recurrence expression. 

7. (Original) The machine-implementable method of claim 5 wherein said 
mathematical recurrence expression is a backward mathematical recurrence expression. 

8. (Original) The machine-implementable method of claim 7 further adapted to 
implement said backward mathematical recurrence expression by comprising further 
steps for: 

e) equating, via said machine-processing unit, a value of a highest ordered 
coefficient of said second polynomial c(x) to a value of an identified highest ordered 
coefficient of said first polynomial p(x) by generating a plurality of machine processing 
unit signals to determine: C a -i=P n ; and 

f) computing, via a machine processing unit, a value for each lower ordered 
coefficient of said second polynomial c(x) by generating a plurality of machine 
processing unit signals to determine:C k -l=Xj • C k +P k where k = (n-1) to 1. 

9. (Original) The machine-implementable method of claim 8 wherein said 
predetermined xi is selected from a set of predetermined values of x*, 

10. (Original) The machine-implementable method of claim 9 wherein said 
predetermined x* is a closest member of said set to said identified x* 
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1 1 . (Original) The machine-implementable method of claim 1 0 wherein said step 
of determining a value of said second polynomial c(x) is computed by using Homer's 
Rule. 

12. (Original) The machine-implementable method of claim 1 1 for determining a 
value of a denominator polynomial q(x) having identified ordered denominator 
coefficients, said denominator polynomial q(x) being expressible as: q(x) « 2(Qj • x 1 ) 
where j=0 to rn, comprising further steps of: 

g) adapting said input data to further include a value for each identified ordered 
denominator coefficient of said denominator polynomial q(x), a value of a predetennined 
q(xO correspondingly paired with said predetermined x t ; and 

h) detennining, via said machine processing unit, a value of another polynomial 
d(x) having ordered denominator coefficients, said another polynomial d(x) being 
expressible as: d(x) = L(Dk • x k ) where k = 0 to (m-1) so that said denominator 
polynomial q(x) is expressible as: q(x) - q(xj)+{x- Xi } * d(x), and a value for the said 
denominator polynomial is resolved. 

13. (Original) The machine-implementable method of claim 12 wherein the first 
polynomial p(x) is a numerator polynomial p(x), and p(x)-p(xO is computed, and p(x>) is 
not read. 
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14. (Original) The machine-implementable method of claim 13 for determining a 

value of a rational function r(x) being expressible as a quotient of said numerator 

polynomial p(x) and said denominator polynomial q(x) expressed as r(x) ta p(x) / q(x), 

comprising further steps of: 

i) adapting said input data to further including a value of a predetermined r(xj) 
correspondingly paired with said predetermined x*; and 

j) constructing, via said machine processing unit, a value of said rational function 
r(x) by generating a plurality of machine processing unit signals to determine: 

r(x) = r(xO • (1 - (q (x)-q (x ; ))/q(x))) + (p(x) - p(xO)/q(x) . 

15. (Original) The machine-implementable method of claim 14 wherein said 
rational function r(x) is an approximation to a Bessel function. 

16. (Original) The machine-implementable method of claim 14 wherein said 
rational function r(x) is an approximation to an error function (ERF). 

17. (Original) The machine-implementable method of claim 14 wherein said 
rational function r(x) is an approximation to a complementary error function (ERFC). 

18. (Original) The machine-implementable method of claim 14 wherein said 
rational function r(x) is an approximation to a log gamma function (LGAMMA). 
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19. (Original) The machine-implementable method of claim 1 1 or 14 wherein all 

values are machine representations of some numerical value, said machine processing 

unit is a computer processing unit, each machine representation is a binary representation 

operable with said computer processing unit, and machine-readable medium is a 

computer-readable medium. 

Claims 20-22 (Cancelled) 

23. (Currently Amended) A machine for computing a property of a 
mathematically modelled physical system, the machine configured to perform the steps 
comprising: 

a) reading, via a machine processing unit, input data including a value for each 
identified ordered coefficient of a first polynomial p(x) representing said property, said 
polynomial p(x) being expressed as p(x) = 2(Pj • x>) where j = 0 to n, a value of a quantity 
x, a value of a predetermined xj, and a value of a predetermined p(x*) correspondingly 
paired with said predetermined x»; 

b) building, via said machine processing unit, a value of a second polynomial c(x) 
having ordered coefficients, said second polynomial c(x) being expressible as: c(x) - 
S(C k • x 14 ) where k - 0 to (n-1) so that said first polynomial p(x) is expressible as: p(x) = 
p(x*)+ {x-xj} • c(x), comprising the steps of: 

i) determining, via said machine processing unit, a value for each ordered 
coefficient of said second polynomial c(x) by generating a plurality of machine 
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processing unit signals to determine each said ordered coefficient of said second 

polynomial c(x) from: C*=2(fVi+j> where j= 0 to (n-l-k); 

ii) detennining, via said machine processing unit, a value of said second 

polynomial c(x) by generating a plurality of machine processing unit signals to 

determine: c(x) = E(Qc • x k ) where k ■* 0 to (n-1); 

c) constructing, via said machine processing unit, a value of said first polynomial 
p(x) by generating a plurality of machine processing unit signals to determine: 
p(x)=p(xi)+{x-Xi} . C (x); and 

d) outputting, via said machine-processing unit, said value of the first polynomial 
p(x) representing said property of the mathematically modelled physical system. 

24. (Original) The machine of claim 23 wherein a difference between x and x\ is 
sufficiently small to achieve a desired accuracy of a final computation of said machine 
representation of a numerical value of said first polynomial p(x). 

25. (Original) The machine of claim 24 wherein said means for reading said 
input data comprises means for reading, via said machine processing unit, said input data 
from a machine-readable medium. 

26. (Original) The machine of claim 25 wherein said ordered coefficients of said 
second polynomial c(x) are computed from a mathematical expression derivable from: 
C fc = I(P ( icH + j) • di) where j= 0 to (n-l-k). 
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27. (Original) The machine of claim 26 wherein said mathematical expression is 
a mathematical recurrence expression. 

28. (Original) The machine of claim 27 wherein said mathematical recurrence 
expression is a forward mathematical recurrence expression. 

29. (Original) The machine of claim 27 wherein said mathematical recurrence 
expression is a backward mathematical recurrence expression. 

30. (Original) The machine of claim 29 further adapted to implement said 
backward mathematical recurrence expression by further comprising: 

e) means for equating, via said machine processing unit, a value of a highest 
ordered coefficient of said second polynomial c(x) to a value of an identified highest 
ordered coefficient of said first polynomial p(x) by generating a plurality of machine 
processing unit signals to determine: C n -i - P ft ; and 

f) means for computing, via said machine processing unit, a value for each lower 
ordered coefficient of said second polynomial c(x) by generating aplurality of machine 
processing unit signals to determine: C k+] = xi • C k +P k where k = (n-1) to 1. 

3 1 . (Original) The machine of claim 30 wherein said predetermined x* is selected 
from a set of predetermined values of Xj, 
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32. (Original) The machine of claim 30 wherein said predetermined Xi is a 
closest member of said set to said identified x. 

33. (Original) The machine of claim 32 wherein the determining means for 
detennining a value of said second polynomial c(x) is computed by using Homer's Rule. 

34. (Original) The machine of claim 33 for determining a value of a denominator 
polynomial q(x) having identified ordered denominator coefficients, said denominator 
polynomial q(x) being expressible as: q(x) = £(Qj • x 1 ) where j = 0 to m, comprising 
further steps of: 

g) adapting said input data to further include a value for each identified ordered 
denominator coefficient of said denominator polynomial q(x), and a value of a 
predetermined q(xj) correspondingly paired with said predetermined x<; and 

h) determining, via said machine processing unit, a value of another polynomial 
d(x) having ordered denominator coefficients, said another polynomial d(x) being 
expressible as: d(x) «= Z(D k • x*) where k 0 to (m-1) so that said denominator 
polynomial q(x) is expressible as: q(x)=q(xj)+{x- xj } • d(x), and a value for the said 
denominator polynomial is resolved, 

35. (Original) The machine of claim 34 wherein the first polynomial p(x) is a 
numerator polynomial p(x), and p(x)-p(x*) is computed, and p(xj) is not read. 
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36. (Original) The machine of claim 35 for determining a value of a rational 
function r(x) being expressible as a quotient of said numerator polynomial p(x) and said 
denominator polynomial q(x) expressed as r(x) = p(x) / q(x), comprising further steps of: 

i) adapting said input data to further including a value of a predetermined r(Xj) 
correspondingly paired with said predetermined x*; and 

j) constructing, via said machine processing unit, a value of said rational function 
r(x) by generating a plurality of machine processing unit signals to determine: 

Hx) - *frO • (1 - (q (x)-q (Xi))/q(x))) + (p(x) - p(xO)/q(x) - 

37. (Original) The machine of claim 36 wherein said rational function is an 
approximation to a Bessel function. 

38. (Original) The machine of claim 36 wherein said rational function is an 
approximation to an error function (ERF). 

39. (Original) The machine of claim 36 wherein said rational function is an 
approximation to a complementary error function (ERFC). 

40. (Original) The machine of claim 36 wherein said rational function is an 
approximation to a log gamma function (LGAMMA). 
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41 . (Original) The machine of claim 33 or 36 wherein all values are machine 
representations of some numerical value, said machine processing unit is a computer 
processing unit, each machine representation is a binary representation operable with said 
computer processing unit, and said machine-readable medium is a computer-readable 
medium. 

42. (Original) A machine having a computer-readable program product having 
computer executable instructions for instructing a computer to embody the machine of 
claim 41. 

43. (Original) A machine having a computer-readable mathematical software 
routine library including computer executable instructions for instructing a computer to 
embody the machine of claim 41 . 

44. (Original) A machine having the computer-readable mathematical software 
routine library of claim 43 wherein said library is operatively associated with a software 
programming language. 
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